Fusion InPipe, an integrative pipeline for gene fusion detection from RNA-seq data in acute pediatric leukemia.

RNA sequencing (RNA-seq) is a reliable tool for detecting gene fusions in acute leukemia. Multiple bioinformatics pipelines have been developed to analyze RNA-seq data, but an agreed gold standard has not been established. This study aimed to compare the applicability of 5 fusion calling pipelines (Arriba, deFuse, CICERO, FusionCatcher, and STAR-Fusion), as well as to define and develop an integrative bioinformatics pipeline (Fusion InPipe) to detect clinically relevant gene fusions in acute pediatric leukemia. We analyzed RNA-seq data by each pipeline individually and by Fusion InPipe. Each algorithm individually called most of the fusions with similar sensitivity and precision. However, not all rearrangements were called, suggesting that choosing a single pipeline might cause missing important fusions. To improve this, we integrated the results of the five algorithms in just one pipeline, Fusion InPipe, comparing the output from the agreement of 5/5, 4/5, and 3/5 algorithms. The maximum sensitivity was achieved with the agreement of 3/5 algorithms, with a global sensitivity of 95%, achieving a 100% in patients' data. Furthermore, we showed the necessity of filtering steps to reduce the false positive detection rate. Here, we demonstrate that Fusion InPipe is an excellent tool for fusion detection in pediatric acute leukemia with the best performance when selecting those fusions called by at least 3/5 pipelines.

Germline mutations in WNK2 could be associated with serrated polyposis syndrome.

BACKGROUND: Patients with serrated polyposis syndrome (SPS) have multiple and/or large serrated colonic polyps and higher risk for colorectal cancer. SPS inherited genetic basis is mostly unknown. We aimed to identify new germline predisposition factors for SPS by functionally evaluating a candidate gene and replicating it in additional SPS cohorts. METHODS: After a previous whole-exome sequencing in 39 SPS patients from 16 families (discovery cohort), we sequenced specific genes in an independent validation cohort of 211 unrelated SPS cases. Additional external replication was also available in 297 SPS cases. The WNK2 gene was disrupted in HT-29 cells by gene editing, and WNK2 variants were transfected using a lentiviral delivery system. Cells were analysed by immunoblots, real-time PCR and functional assays monitoring the mitogen-activated protein kinase (MAPK) pathway, cell cycle progression, survival and adhesion. RESULTS: We identified 2 rare germline variants in the WNK2 gene in the discovery cohort, 3 additional variants in the validation cohort and 10 other variants in the external cohorts. Variants c.2105C>T (p.Pro702Leu), c.4820C>T (p.Ala1607Val) and c.6157G>A (p.Val2053Ile) were functionally characterised, displaying higher levels of phospho-PAK1/2, phospho-ERK1/2, CCND1, clonogenic capacity and MMP2. CONCLUSION: After whole-exome sequencing in SPS cases with familial aggregation and replication of results in additional cohorts, we identified rare germline variants in the WNK2 gene. Functional studies suggested germline WNK2 variants affect protein function in the context of the MAPK pathway, a molecular hallmark in this disease.

Long-Term Strenuous Exercise Promotes Vascular Injury by Selectively Damaging the Tunica Media: Experimental Evidence.

Moderate exercise has well-founded benefits in cardiovascular health. However, increasing, yet controversial, evidence suggests that extremely trained athletes may not be protected from cardiovascular events as much as moderately trained individuals. In our rodent model, intensive but not moderate training promoted aorta and carotid stiffening and elastic lamina ruptures, tunica media thickening of intramyocardial arteries, and an imbalance between vasoconstrictor and relaxation agents. An up-regulation of angiotensin-converter enzyme, miR-212, miR-132, and miR-146b might account for this deleterious remodeling. Most changes remained after a 4-week detraining. In conclusion, our results suggest that intensive training blunts the benefits of moderate exercise.

Molecular Modelling Hurdle in the Next-Generation Sequencing Era.

There are challenges in the genetic diagnosis of rare diseases, and pursuing an optimal strategy to identify the cause of the disease is one of the main objectives of any clinical genomics unit. A range of techniques are currently used to characterize the genomic variability within the human genome to detect causative variants of specific disorders. With the introduction of next-generation sequencing (NGS) in the clinical setting, geneticists can study single-nucleotide variants (SNVs) throughout the entire exome/genome. In turn, the number of variants to be evaluated per patient has increased significantly, and more information has to be processed and analyzed to determine a proper diagnosis. Roughly 50% of patients with a Mendelian genetic disorder are diagnosed using NGS, but a fair number of patients still suffer a diagnostic odyssey. Due to the inherent diversity of the human population, as more exomes or genomes are sequenced, variants of uncertain significance (VUSs) will increase exponentially. Thus, assigning relevance to a VUS (non-synonymous as well as synonymous) in an undiagnosed patient becomes crucial to assess the proper diagnosis. Multiple algorithms have been used to predict how a specific mutation might affect the protein's function, but they are far from accurate enough to be conclusive. In this work, we highlight the difficulties of genomic variability determined by NGS that have arisen in diagnosing rare genetic diseases, and how molecular modelling has to be a key component to elucidate the relevance of a specific mutation in the protein's loss of function or malfunction. We suggest that the creation of a multi-omics data model should improve the classification of pathogenicity for a significant amount of the detected genomic variability. Moreover, we argue how it should be incorporated systematically in the process of variant evaluation to be useful in the clinical setting and the diagnostic pipeline.

Muscular carnosine is a marker for cardiorespiratory fitness and cardiometabolic risk factors in men with type 1 diabetes.

PURPOSE: Muscle is an essential organ for glucose metabolism and can be influenced by metabolic disorders and physical activity. Elevated muscle carnosine levels have been associated with insulin resistance and cardiometabolic risk factors. Little is known about muscle carnosine in type 1 diabetes (T1D) and how it is influenced by physical activity. The aim of this study was to characterize muscle carnosine in vivo by proton magnetic resonance spectroscopy (1H MRS) and evaluate the relationship with physical activity, clinical characteristics and lipoprotein subfractions. METHODS: 16 men with T1D (10 athletes/6 sedentary) and 14 controls without diabetes (9/5) were included. Body composition by DXA, cardiorespiratory capacity (VO2peak) and serum lipoprotein profile by proton nuclear magnetic resonance (1H NMR) were obtained. Muscle carnosine scaled to water (carnosineW) and to creatine (carnosineCR), creatine and intramyocellular lipids (IMCL) were quantified in vivo using 1H MRS in a 3T MR scanner in soleus muscle. RESULTS: Subjects with T1D presented higher carnosine CR levels compared to controls. T1D patients with a lower VO2peak presented higher carnosineCR levels compared to sedentary controls, but both T1D and control groups presented similar levels of carnosineCR at high VO2peak levels. CarnosineW followed the same trend. Integrated correlation networks in T1D demonstrated that carnosineW and carnosineCR were associated with cardiometabolic risk factors including total and abdominal fat, pro-atherogenic lipoproteins (very low-density lipoprotein subfractions), low VO2peak, and IMCL. CONCLUSIONS: Elevated muscle carnosine levels in persons with T1D and their effect on atherogenic lipoproteins can be modulated by physical activity.

The Increasing Impact of Translational Research in the Molecular Diagnostics of Neuromuscular Diseases.

The diagnosis of neuromuscular diseases (NMDs) has been progressively evolving from the grouping of clinical symptoms and signs towards the molecular definition. Optimal clinical, biochemical, electrophysiological, electrophysiological, and histopathological characterization is very helpful to achieve molecular diagnosis, which is essential for establishing prognosis, treatment and genetic counselling. Currently, the genetic approach includes both the gene-targeted analysis in specific clinically recognizable diseases, as well as genomic analysis based on next-generation sequencing, analyzing either the clinical exome/genome or the whole exome or genome. However, as of today, there are still many patients in whom the causative genetic variant cannot be definitely established and variants of uncertain significance are often found. In this review, we address these drawbacks by incorporating two additional biological omics approaches into the molecular diagnostic process of NMDs. First, functional genomics by introducing experimental cell and molecular biology to analyze and validate the variant for its biological effect in an in-house translational diagnostic program, and second, incorporating a multi-omics approach including RNA-seq, metabolomics, and proteomics in the molecular diagnosis of neuromuscular disease. Both translational diagnostics programs and omics are being implemented as part of the diagnostic process in academic centers and referral hospitals and, therefore, an increase in the proportion of neuromuscular patients with a molecular diagnosis is expected. This improvement in the process and diagnostic performance of patients will allow solving aspects of their health problems in a precise way and will allow them and their families to take a step forward in their lives.

Maternal proteomic profiling reveals alterations in lipid metabolism in late-onset fetal growth restriction.

Fetal growth restriction defined as the failure to achieve the fetal genetic growth potential is a major cause of perinatal morbidity and mortality. The role of maternal adaptations to placental insufficiency in this disorder is still not fully understood. We aimed to investigate the biological processes and protein-protein interactions involved in late-onset fetal growth restriction in particular. We applied 2D nano LC-MS/MS proteomics analysis on maternal blood samples collected at the time of delivery from 5 singleton pregnancies with late-onset fetal growth restriction and 5 uncomplicated pregnancies. Data were analyzed using R package "limma" and Ingenuity Pathway Analysis. 25 proteins showed significant changes in their relative abundance in late-onset fetal growth restriction (p value < 0.05). Direct protein-protein interactions network demonstrated that Neurogenic locus notch homolog protein 1 (NOTCH1) was the most significant putative upstream regulator of the observed profile. Gene ontology analysis of these proteins revealed the involvement of 14 canonical pathways. The most significant biological processes were efflux of cholesterol, efflux of phospholipids, adhesion of blood cells, fatty acid metabolism and dyslipidemia. Future studies are warranted to validate the potential role of the detected altered proteins as potential therapeutic targets in the late-onset form of fetal growth restriction.

Helpful Criteria When Implementing NGS Panels in Childhood Lymphoblastic Leukemia.

The development of Next-Generation Sequencing (NGS) has provided useful diagnostic, prognostic, and therapeutic strategies for individualized management of B-cell precursor acute lymphoblastic leukemia (BCP-ALL) patients. Consequently, NGS is rapidly being established in clinical practice. However, the technology's complexity, bioinformatics analysis, and the different available options difficult a broad consensus between different laboratories in its daily routine introduction. This collaborative study among Spanish centers was aimed to assess the feasibility, pros, and cons of our customized panel and other commercial alternatives of NGS-targeted approaches. The custom panel was tested in three different sequencing centers. We used the same samples to assess other commercial panels (OncomineTM Childhood Cancer Research Assay; Archer®FusionPlex® ALL, and Human Comprehensive Cancer Panel GeneRead Panel v2®). Overall, the panels showed a good performance in different centers and platforms, but each NGS approach presented some issues, as well as pros and cons. Moreover, a previous consensus on the analysis and reporting following international guidelines would be preferable to improve the concordance in results among centers. Our study shows the challenges posed by NGS methodology and the need to consider several aspects of the chosen NGS-targeted approach and reach a consensus before implementing it in daily practice.

Coenzyme Q10 Treatment Monitoring in Different Human Biological Samples.

Coenzyme Q10 (CoQ) treatment monitoring is a matter of debate since CoQ distribution from plasma to blood cells and tissues is not fully understood. We aimed to analyze the CoQ levels in a wide set of human biological samples (plasma, blood mononuclear cells (BMCs), platelets, urinary cells, and skeletal muscle) from a group of 11 healthy male runners before and after CoQ supplementation. The CoQ content in the different samples was analyzed by HPLC coupled to electrochemical detection. No significant differences were observed in the CoQ levels measured in the BMCs, platelets, and urine after the one-month treatment period. Plasma CoQ (expressed in absolute values and values relative to total cholesterol) significantly increased after CoQ supplementation (p = 0.003 in both cases), and the increase in CoQ in muscle approached significance (p = 0.074). CoQ levels were increased in the plasma of all supplemented subjects, and muscle CoQ levels were increased in 8 out of 10 supplemented subjects. In conclusion, the analysis of CoQ in plasma samples seems to be the best surrogate biomarker for CoQ treatment monitoring. Moreover, oral CoQ administration was effective for increasing muscle CoQ concentrations in most subjects.

Integrated mRNA and microRNA transcriptome profiling during differentiation of human nasal polyp epithelium reveals an altered ciliogenesis.

BACKGROUND: Human adult basal stem/progenitor cells (BSCs) obtained from chronic rhinosinusitis with nasal polyps (CRSwNP) when differentiated in an air-liquid interface (ALI) usually provide a pseudostratified airway epithelium with similar abnormalities than original in vivo phenotype. However, the intrinsic mechanisms regulating this complex process are not well defined and their understanding could offer potential new therapies for CRSwNP (incurable disease). METHODS: We performed a transcriptome-wide analysis during in vitro mucociliary differentiation of human adult BSCs from CRSwNP, compared to those isolated from control nasal mucosa (control-NM), in order to identify which key mRNA and microRNAs are regulating this complex process in pathological and healthy conditions. RESULTS: A number of genes, miRs, biological processes, and pathways were identified during mucociliary differentiation of both CRSwNP and control-NM epithelia, and notably, we have demonstrated for the first time that genetic transcriptional program responsible of ciliogenesis and cilia function is significantly impaired in CRSwNP epithelium, presumably produced by an altered expression of microRNAs, particularly of those miRs belonging to mir-34 and mi-449 families. CONCLUSIONS: This study provides for the first time a novel insight into the molecular basis of sinonasal mucociliary differentiation, demonstrating that transcriptome related to ciliogenesis and cilia function is significantly impaired during differentiation of CRSwNP epithelium due to an altered expression of microRNAs.

Comprehensive Analysis of GABAA-A1R Developmental Alterations in Rett Syndrome: Setting the Focus for Therapeutic Targets in the Time Frame of the Disease.

Rett syndrome, a serious neurodevelopmental disorder, has been associated with an altered expression of different synaptic-related proteins and aberrant glutamatergic and γ-aminobutyric acid (GABA)ergic neurotransmission. Despite its severity, it lacks a therapeutic option. Through this work we aimed to define the relationship between MeCP2 and GABAA.-A1 receptor expression, emphasizing the time dependence of such relationship. For this, we analyzed the expression of the ionotropic receptor subunit in different MeCP2 gene-dosage and developmental conditions, in cells lines, and in primary cultured neurons, as well as in different developmental stages of a Rett mouse model. Further, RNAseq and systems biology analysis was performed from post-mortem brain biopsies of Rett patients. We observed that the modulation of the MeCP2 expression in cellular models (both Neuro2a (N2A) cells and primary neuronal cultures) revealed a MeCP2 positive effect on the GABAA.-A1 receptor subunit expression, which did not occur in other proteins such as KCC2 (Potassium-chloride channel, member 5). In the Mecp2+/- mouse brain, both the KCC2 and GABA subunits expression were developmentally regulated, with a decreased expression during the pre-symptomatic stage, while the expression was variable in the adult symptomatic mice. Finally, the expression of the gamma-aminobutyric acid (GABA) receptor-related synaptic proteins from the postmortem brain biopsies of two Rett patients was evaluated, specifically revealing the GABA A1R subunit overexpression. The identification of the molecular changes along with the Rett syndrome prodromic stages strongly endorses the importance of time frame when addressing this disease, supporting the need for a neurotransmission-targeted early therapeutic intervention.

Effect of blood contamination of cerebrospinal fluid on amino acids, biogenic amines, pterins and vitamins.

BACKGROUND: Cerebrospinal fluid (CSF) metabolomic investigations are a powerful tool for studying neurometabolic diseases. We aimed to assess the effect of CSF contamination with blood on the concentrations of selected biomarkers. METHODS: CSF samples were spiked in duplicate with increasing volumes of whole blood under two conditions: (A) pooled CSF spiked with fresh whole blood and frozen to cause red blood cell (RBC) lysis; (B) pooled CSF spiked with fresh blood and centrifuged (the supernatant with no RBCs was frozen until the moment of analysis). CSF concentrations of amino acids, biogenic amines, pterins, and vitamins were analysed by HPLC coupled with tandem mass spectrometry, electrochemical and fluorescence detection. RESULTS: Aspartate, glutamate, taurine, ornithine, glycine, citrulline, pyridoxal 5´-phosphate, 5-methyltetrahydrofolate, and thiamine showed higher values when RBCs were lysed when compared with those of CSF with no RBC, while arginine, 5-hydroxyindoleacetic and homovanillic acids showed lower values. When RBCs were removed from CSF, only some amino acids, thiamine and pyridoxal 5´-phosphate showed moderately higher values when compared with the non-spiked CSF sample. CONCLUSIONS: CSF-targeted metabolomic analysis is feasible even when substantial RBC contamination of CSF has occurred since CSF centrifugation to remove RBC prior to freezing eliminated most of the interferences observed.

Plasma coenzyme Q10 status is impaired in selected genetic conditions.

Identifying diseases displaying chronic low plasma Coenzyme Q10 (CoQ) values may be important to prevent possible cardiovascular dysfunction. The aim of this study was to retrospectively evaluate plasma CoQ concentrations in a large cohort of pediatric and young adult patients. We evaluated plasma CoQ values in 597 individuals (age range 1 month to 43 years, average 11 years), studied during the period 2005-2016. Patients were classified into 6 different groups: control group of healthy participants, phenylketonuric patients (PKU), patients with mucopolysaccharidoses (MPS), patients with other inborn errors of metabolism (IEM), patients with neurogenetic diseases, and individuals with neurological diseases with no genetic diagnosis. Plasma total CoQ was measured by reverse-phase high-performance liquid chromatography with electrochemical detection and ultraviolet detection at 275 nm. ANOVA with Bonferroni correction showed that plasma CoQ values were significantly lower in the PKU and MPS groups than in controls and neurological patients. The IEM group showed intermediate values that were not significantly different from those of the controls. In PKU patients, the Chi-Square test showed a significant association between having low plasma CoQ values and being classic PKU patients. The percentage of neurogenetic and other neurological patients with low CoQ values was low (below 8%). In conclusión, plasma CoQ monitoring in selected groups of patients with different IEM (especially in PKU and MPS patients, but also in IEM under protein-restricted diets) seems advisable to prevent the possibility of a chronic blood CoQ suboptimal status in such groups of patients.

Association of T and NK Cell Phenotype With the Diagnosis of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS).

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a pathological condition characterized by incapacitating fatigue and a combination of neurologic, immunologic, and endocrine symptoms. At present its diagnosis is based exclusively on clinical criteria. Several studies have described altered immunologic profiles; therefore, we proposed to further examine the more significant differences, particularly T and NK cell subpopulations that could be conditioned by viral infections, to discern their utility in improving the diagnosis and characterization of the patients. The study included 76 patients that fulfilled the revised Canadian Consensus Criteria (CCC 2010) for ME/CFS and 73 healthy controls, matched for age and gender. Immunophenotyping of different T cell and natural killer cell subpopulations in peripheral blood was determined by flow cytometry. ME/CFS patients showed significantly lower values of T regulatory cells (CD4+CD25++(high)FOXP3+) and higher NKT-like cells (CD3+CD16+/-CD56+) than the healthy individuals. Regarding NK phenotypes, NKG2C was significantly lower and NKCD69 and NKCD56 bright were significantly higher in the patients group. A classification model was generated using the more relevant cell phenotype differences (NKG2C and T regulatory cells) that was able to classify the individuals as ME/CFS patients or healthy in a 70% of cases. The observed differences in some of the subpopulations of T and NK cells between patients and healthy controls could define a distinct immunological profile that can help in the diagnostic process of ME/CFS patients, contribute to the recognition of the disease and to the search of more specific treatments. However, more studies are needed to corroborate these findings and to contribute to establish a consensus in diagnosis.

Characterization of Inflammatory Response in Acute-on-Chronic Liver Failure and Relationship with Prognosis.

ACLF is characterized by a systemic inflammatory response, but the cytokines involved in this process have not been well studied. The aim of this study was to characterize the systemic inflammatory response in patients with cirrhosis and ACLF and its relationship with prognosis. Fifty-five patients with cirrhosis, 26 with ACLF, were studied prospectively. Systemic inflammatory response was analyzed by measuring a large array of plasma cytokines by using a multiplex kit. A principal component analysis show noticeable differences between ACLF and decompensated cirrhosis without ACLF. Patients with ACLF had significant abnormal levels of 12 cytokines compared to those without ACLF, including: VCAM-1, VEGF-A, Fractalkine, MIP-1α, Eotaxin, IP-10, RANTES, GM-CSF, IL-1ß, IL-2, ICAM-1, and MCP-1. Cytokines showing the most marked relationship with ACLF were VCAM-1 and VEGF-A (AUCROC 0.77; p = 0.001). There was a significant relationship between some of inflammatory mediators and 3-month mortality, particularly VCAM-1, ICAM-1, and GM-CSF (AUCROC>0.7; p < 0.05). Functional Enrichment Analysis showed that inflammatory markers differentially expressed in ACLF patients were enriched in leukocyte migration, particularly monocytes and macrophages, and chemotaxis pathways. In conclusion, ACLF is characterized by a marked inflammatory reaction with activation of mediators of adhesion and migration of leukocytes. The intensity of the inflammatory reaction correlates with prognosis.

Targeting IκB proteins for HIV latency activation: the role of individual IκB and NF-κB proteins.

Latently infected cell reservoirs represent the main barrier to HIV eradication. Combination antiretroviral therapy (cART) effectively blocks viral replication but cannot purge latent provirus. One approach to HIV eradication could include cART to block new infections plus an agent to activate latent provirus. NF-κB activation induces HIV expression, ending latency. Before activation, IκB proteins sequester NF-κB dimers in the cytoplasm. Three canonical IκBs, IκBα, IκBß, and IκBε, exist, but the IκB proteins' role in HIV activation regulation is not fully understood. We studied the effects on HIV activation of targeting IκBs by single and pairwise small interfering RNA (siRNA) knockdown. After determining the relative abundance of the IκBs, the relative abundance of NF-κB subunits held by the IκBs, and the kinetics of IκB degradation and resynthesis following knockdown, we studied HIV activation by IκB knockdown, in comparison with those of known HIV activators, tumor necrosis factor alpha (TNF-α), tetradecanoyl phorbol acetate (TPA), and trichostatin A (TSA), in U1 monocytic and J-Lat 10.6 lymphocytic latently infected cells. We found that IκBα knockdown activated HIV in both U1 and J-Lat 10.6 cells, IκBß knockdown did not activate HIV, and, surprisingly, IκBε knockdown produced the most HIV activation, comparable to TSA activation. Our data show that HIV reactivation can be triggered by targeting two different IκB proteins and that IκBε may be an effective target for HIV latency reactivation in T-cell and macrophage lineages. IκBε knockdown may offer attractive therapeutic advantages for HIV activation because it is not essential for mammalian growth and development and because new siRNA delivery strategies may target siRNAs to HIV latently infected cells.

Cell line-dependent variability in HIV activation employing DNMT inhibitors.

Long-lived reservoirs of Human Immunodeficiency Virus (HIV) latently infected cells present the main barrier to a cure for HIV infection. Much interest has focused on identifying strategies to activate HIV, which would be used together with antiretrovirals to attack reservoirs. Several HIV activating agents, including Tumor Necrosis Factor alpha (TNFα) and other agents that activate via NF-kB are not fully effective in all latent infection models due to epigenetic restrictions, such as DNA methylation and the state of histone acetylation. DNA methyltransferases (DNMT) inhibitors like 5-aza-2'deoxycytidine (Aza-CdR) and histone deacetylase (HDAC) inhibitors like Trichostatin A (TSA) have been proposed as agents to enhance reactivation and have shown activity in model systems. However, it is not clear how the activities of DNMT and HDAC inhibitors range across different latently infected cell lines, potential models for the many different latently infected cells within an HIV patient. We determined HIV activation following treatment with TNFα, TSA and Aza-CdR across a range of well known latently infected cell lines. We assessed the activity of these compounds in four different Jurkat T cell-derived J-Lat cell lines (6.3, 8.4, 9.2 and 10.6), which have a latent HIV provirus in which GFP replaces Nef coding sequence, and ACH-2 and J1.1 (T cell-derived), and U1 (promonocyte-derived) cell lines with full-length provirus. We found that Aza-CdR plus TNFα activated HIV at least twice as well as TNFα alone for almost all J-Lat cells, as previously described, but not for J-Lat 10.6, in which TNFα plus Aza-CdR moderately decreased activation compared to TNFα alone. Surprisingly, a much greater reduction of TNFα-stimulated activation with Aza-CdR was detected for ACH-2, J1.1 and U1 cells. Reaching the highest reduction in U1 cells with a 75% reduction. Interestingly, Aza-CdR not only decreased TNFα induction of HIV expression in certain cell lines, but also decreased activation by TSA. Since DNMT inhibitors reduce the activity of provirus activators in some HIV latently infected cell lines the use of epigenetic modifying agents may need to be carefully optimized if they are to find clinical utility in therapies aimed at attacking latent HIV reservoirs.

Analysis of recombinant protein toxicity in E. coli through a phage lambda-based genetic screening system.

The aspartic protease from the human immunodeficiency virus type 1 (HIV-1) is highly toxic to E. coli, thus impairing its yield in production processes. Proteolytic cleavage of essential cellular proteins is probably a major contributor to the bacteriocidal effect but this has not been proven. Through an adapted high-throughput lambda-based screening system, we have analyzed a set of HIV-1 protease mutants with distinguishable catalytic properties and we show that inactive enzymes are as toxic to E. coli cells as the wild-type enzyme. Together with additional data from directed molecular evolution approaches, these results indicate that the toxicity of the viral protease is not linked to its proteolytic activity. Our study also reveals that the lambda-based screening system is a robust new tool for the genetic analysis of highly toxic recombinant products in E. coli.

HIV-1 protease catalytic efficiency effects caused by random single amino acid substitutions.

Protein evolution has occurred by successive fixation of individual mutations. The probability of fixation depends on the fitness of the mutation, and the arising variant can be deleterious, neutral, or beneficial. Despite its relevance, only few studies have estimated the distribution of fitness effects caused by random single mutations on protein function. The human immunodeficiency virus type 1 (HIV-1) protease was chosen as a model protein to quantify protein's tolerability to random single mutations. After determining the enzymatic activity of 107 single random mutants, we found that 86% of single mutations were deleterious for the enzyme catalytic efficiency and 54% lethal. Only 2% of the mutations significantly increased the catalytic efficiency of the enzyme. These data demonstrate the vulnerability of HIV-1 protease to single random mutations. When a second random mutagenesis library was constructed from an HIV-1 protease carrying a highly deleterious single mutation (D30N), a higher proportion of mutations with neutral or beneficial effect were found, 26% and 9%, respectively. Importantly, antagonist epistasis was observed between deleterious mutations. In particular, the mutation N88D, lethal for the wild-type protease, restored the wild-type catalytic efficiency when combined with the highly deleterious mutation D30N. The low tolerability to single random substitutions shown here for the wild-type HIV-1 protease contrasts with its in vivo ability to generate an adaptive variation. Thus, the antagonist epistasis between deleterious or lethal mutations may be responsible for increasing the protein mutational robustness and evolvability.